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Abstract 

We consider M = 2 quiver Chern-Simons theories described by 
brane tilings, whose moduli spaces are toric Calabi-Yau 4-folds. Sim- 
ple prescriptions to obtain toric data of the moduli space and a cor- 
responding brane crystal from a brane tiling are proposed. 
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1 Introduction 



Recently, three-dimensional supersymmetric Chern-Simons theories have attracted 
great interest as theories for multiple M2-branes in various backgrounds. This 
was triggered by the proposal of M = 8 interacting Chern-Simons theories by 
Bagger, Lambert [U [21 [3], and GusstavsonJU [5]. Their model (BLG model) is 
based on Lie 3-algebra, and the action includes structure constants, which sat- 
isfy the so-called fundamental identity. This model, however, has not succeeded 
in describing an arbitrary number of M2-branes in uncompactified flat back- 
ground, due to the fact that the fundamental identity is very restrictive and it 
admits the only one non-trivial finite-dimensional algebra with a positive definite 
metric[6l [7]. The resulting theory is conjectured to describe two M2-branes on a 
certain orbifoldjHJ E]. 

Aharony et al. proposed alternative model in [TO], based on the recent 
progress in M = 4 Chern-Simons theories [TT1 [T2] . Their model (ABJM model) is 
U(N) x U(N) Chern-Simons gauge theory at level (k, —k) with bi-fundamental 
matter fields. The model describes N M2-branes in the C 4 /Z^ orbifold back- 
ground. Although only M = 6 supersymmetry is manifest in the model, the 
supersymmetry is expected somehow to be enhanced to Af = 8 when k = 1,2. 

After the proposal of the ABJM model, some generalizations have been stud- 
ied. Orbifolds of the ABJM model are discussed in [131 US] • In [H], a certain 
class of M = 3 quiver Chern-Simons theories with non-toric moduli spaces are also 
studied based on the brane construction, and the hyper-Kahler toric structure of 
the moduli spaces is clarified in [16] . The moduli spaces of other superconformal 
Chern-Simons theories are studied in [T71 dU [191 EQl l2T] . 

M = 2 Chern-Simons theories with general quiver structure are studied in [22j, 
and it is shown how the gauge symmetries and D-term conditions are modified 
compared to the case of four- dimensional Af = 1 gauge theories described by the 
same quiver diagrams. It is also found that the moduli spaces of such theories 
generically include a baryonic branch. 

In the case of four- dimensional M = 1 supersymmetric gauge theories, brane 
tilings [231 124"! [25] are convenient tools to establish the relation between gauge 
theories and their moduli spaces, for a class of theories whose moduli spaces are 
toric Calabi-Yau 3-folds. See [26, |27] for review of brane tilings. Brane tilings 
are expected to be convenient for three-dimensional Chern-Simons theory, too. 
M = 2 quiver Chern-Simons theories described by brane tilings are studied in 
[28], and it is shown that the moduli space of the theories are toric Calabi-Yau 
4-folds, and the Hilbert series is computed for some examples. 

In this paper we consider the class of Af = 2 quiver Chern-Simons theories 
described by brane tilings. Our aim is to establish the relation between brane 
tilings and brane crystals. Brane crystals are three-dimensional graphs proposed 
in [291 ESS ED] as diagrams describing three-dimensional superconformal field the- 
ories and the structure of their moduli spaces. We first give a simple prescription 
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to obtain toric data of the moduli space from a tiling, and explain how we can 
construct a crystal describing the same moduli space. 

This paper is organized as follows. In the next section, we explain the relation 
between brane tilings and quiver Chern-Simons theories. In section [3] we review 
how gauge symmetries of Chern-Simons theories are broken due to the existence 
of Chern-Simons terms following [22J. In section HI we define gauge invariant 
operators which parameterize the moduli spaces of Chern-Simons theories. In 
section [5] we give a simple prescription to obtain the toric data of the moduli 
spaces by using tilings. This section has some overlap with [32]. The relation 
between brane tilings and brane crystals are discussed in El The last section is 
devoted to conclusions. 

2 Tilings and Chern-Simons theories 

We consider three-dimensional M = 2 quiver Chern-Simons theories described by 
brane tilings, which are also studied in [28] . 

A brane tiling is a bipartite graph drawn on T 2 . A bipartite graph is a graph 
consisting of vertices of two colors, say, white and black, and all links connect two 
vertices with different colors. Tilings have been used to describe four-dimensional 
Af = 1 quiver gauge theories and the structure of their moduli spaces. The gauge 
group, the matter content, and the superpotential of a gauge theory can be 
read off from the brane tiling for the theory. Namely, faces correspond to U(N) 
factors in the gauge group, and links to bi-fundamental fields. The superpotential 
can be also read off from the tiling in the way we will mention later. These 
correspondences are naturally understood by regarding the tiling as a NS5-D5 
system in type IIB string theory. 

In this paper, we use tilings to describe three-dimensional M = 2 Chern- 
Simons theories. The gauge group, the matter content, and the superpotential 
are read off from the tiling in the same way as the four-dimensional case. These 
rules are naturally understood by regarding the tiling as a D4-NS5 system in type 
IIA theory, rather than the type IIB brane system. By this reason, when we want 
to specify which of three or four dimensional theory a brane tiling describes, we 
call it IIB tiling (for four- dimensional theory), or IIA tiling (for Chern-Simons 
theory). Figure [1] shows an example of IIA tiling for the ABJM model. 

Because we are here interested in the structure of the background spacetime 
probed by M2-branes, we discuss only abelian (N = 1) case. We use indices 
I, J, . . . for links, i, j, . . . for faces, and a,b, . . . for vertices. Let U be the gauge 
group associated with face i, and $/ be the bi-fundamental field associated with 
link I. I £ i means the link / is on the face i. The bi-fundamental field is 
charged under two U(l) factors corresponding to the two faces sharing the link. 
The U(l)i charge Qn of the chiral multiplet $/ is uniquely determined by the 
bipartite graph. When the link / is not a side of the face i Qn = 0. Qn is +1 
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Figure 1: The tiling for the ABJM model at level (k, —k). The arrow represents 
the flow s defined in (J2J). 

(— 1) if the left endpoint of the link / is black (white) when it is seen from the 
face i. 

In order to specify a Chern-Simons theory, we need to fix the levels hi G Z 
for each gauge group as numbers assigned to faces in a IIA tiling. As is pointed 
out in [22], we need to impose the condition 

£** = o, (i) 

i 

to obtain a four- dimensional moduli space. Because of this condition we can 
represent the levels as 

k i = ^2Qns I . (2) 
I 

In order to interpret relations like (j2J) geometrically, we define two kinds of 
flows on the tiling. Let // be a set of numbers assigned to links. A normal flow 
/ is a flow from faces to faces. We define the orientation of the flow to be the 
anti-clockwise direction around black vertices, ((a) in Figure [2]) The other flow 




Figure 2: For a set of numbers // assigned to links we define two flows, (a) is a 
normal flow / and (b) is a tangential flow /*. These two flows are related by the 
7r/2 rotation of arrows. 

associated with // is the tangential flow /*, which describes flow along links from 
black vertices to white ones, ((b) in Figure |2J) We can rewrite the relation in 
terms of the normal flow s or the tangential flow s* as 

{ki} = div s = rot s*. (3) 
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We will later see that this relation has natural interpretation in the context of 
brane realization of Chern-Simons theories. 

For later convenience, we introduce the following normal flows. (See Figure 

El) 

• a : a unit flow along a-cycle on the torus. 

• f3 : a unit flow along j3- cycle on the torus. 

• 7 a : a unit flow around vertex a. The orientation is anti-clockwise (clock- 
wise) around black (white) vertices. 




Figure 3: Examples of cycles a, f3 and 7 a for the ABJM tiling are shown. 

An arbitrary conserved flow can be given as a linear combination of these flows. 
We define an operator Of for a normal flux / with non-negative integral compo- 
nents // by 

O f = W- ( 4 ) 

7 

If / is conserved flow satisfying div / = 0, the operator is also defined for general 
N as a single or multiple trace operator, which is often called mesonic operators. 
Baryonic operators are associated with non-conserved flows. With this notation, 
the superpotential is represented as 

W = £±tt», (5) 

a 

where the signature of the summand is positive (negative) for black (white) ver- 
tices. 

3 Gauge symmetries 

In both three- and four-dimensional cases, the moduli space is defined as the 
coset X/G, where X is the manifold defined by the F-term conditions and G is 
the complexified gauge group. Because IIA and IIB tilings give the same F-term 
conditions, the manifold X is common to two cases. A difference arises in the 
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gauge symmetry G. In this section we review how this difference arises following 
[22]. 

If the tiling has n faces, there are n U(l) factors. Among them, the diagonal 
£7(1) decouples from matter fields, and the effective gauge symmetry is U(l) n ~ l . 
The complexification of this symmetry gives G in the IIB case. In the IIA case, 
however, it is known that U{l) n ~ 1 is broken down to U(l) n ~ 2 due to the existence 
of the Chern-Simons terms. 

Let Ai be the U(l)i gauge field. We define gauge fields 

n n 

a = Y,A, b = Y J KA l . (6) 

1=1 8=1 

Let Cfc (k = 1, . . . , n — 2) be linear comminations of Ai linearly independent of a 
and b. We can rewrite the Chern-Simons terms in the form 

Scs = ^JbAf + S'[b,c k ] (7) 

where f = da and S' does not depend on the diagonal U(l) gauge field a. Because 
the gauge field a does not couple to matter fields, it appears only in the first term 
of (I7j). The action includes a only through /, and we can dualize it by introducing 
Lagrange multiplier r and adding the following term to the action. 

Sr = ~^I dTAf - (8) 
The equation of motion for / is 

dr = -b. (9) 

n 

Let us consider gauge transformation 

Ai -»■ Ai + dOi. (10) 
The relation (jUJ) implies that the dual scalar field r should be transformed under 

<m by 

1 n 

5r = -J2k l 6 l . (11) 

This non-linear gauge transformation of r means that the gauge symmetry is 
always partially broken due to the vev of the scalar field r. As is shown in [22] 
the period of r is 2ir/n, and the parameters for unbroken gauge transformations 
should satisfy 

n 

2nZ3Y,k i 9 i = Y,QiiSi0i> (12) 

i=l I,i 

where we used (J2J to obtain the final expression. 
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4 Gauge invariant operators 



When we analyze the moduli space of a gauge theory, it is convenient to use 
gauge invariant operators as coordinates of the moduli space. In four-dimensional 
gauge theories described by IIB brane tiling, it is known that the moduli space 
is parameterized by three gauge invariant operators 

M a = O a , Mp = O f3 , W = 0~ fa , (13) 

associated with the flows defined in section [2j For general N, these operators 
are defined as single-trace mesonic operators. W is one of the terms in the 
superpotential (jSJ). Due to the F-term conditions, all terms in the superpotential 
have the same vev, and W as an element of the chiral ring does not depend on the 
choice of the vertex a. Because ck, (3, and *y a generate arbitrary conserved flows, 
an arbitrary mesonic operator can be written as a function of these mesonic 
operators, and we can use these three as coordinates in the three-dimensional 
moduli space of four- dimensional gauge theory. 

In the case of Chern-Simons theory, ( TTBl) are again gauge invariant operators, 
and we can use them as coordinates in the moduli space. However, we need 
another gauge invariant operator to parameterize the four- dimensional moduli 
space. Indeed, the restriction (Tl2|) of the gauge transformation parameters admits 
extra gauge invariant operators in addition to the mesonic operators in the four- 
dimensional gauge theory. 

Let us consider an operator O q associated with a flow q, which is not nec- 
essarily conserved. The field $/ associated with link I is transformed under the 
gauge transformation (ITU1) as 

<Z>j^exp(iJ2QnO^$i, (14) 
and the gauge transformation of the operator O q is 

O q ^exp (iJ^Qnqtf^J O q . (15) 

For the operator to be gauge invariant, the components qi of the flow must satisfy 

2nZ 3 J2QuQA- (16) 

If this condition were imposed for arbitrary 8i, solutions would be given by qi = ci 
where c/ is an arbitrary flow satisfying 

Y,QnCi = Vz. (17) 
i 
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This is equivalent to div c = 0, and the normal flow c is conserved. Solutions in 
this form correspond to mesonic operators generated by (TTBl . 

The parameters Q{ are, however, constrained by (fl2i) in the Chern-Simons 
theory. Thus we have an extra solution qj = sj, and a general solution is given 
by 

qj = msj + Cj, (18) 

where m is an arbitrary integer and cj is a conserved flow satisfying (fT7|) . There- 
fore, as the fourth coordinate on the moduli space, we should introduce the 
following baryonic operator associated with s. 

B = 0[s\=H$ a I I . (19) 

An arbitrary gauge invariant operator in the Chern-Simons theory is given as a 
function of the four operators 

M a , M , B, W. (20) 



5 Toric data 

In this section, we give a simple prescription to obtain toric data of the moduli 
space of Chern-Simons theory from the IIA brane tiling for the theory. The same 
subject is also investigated in [32J. 

In general a toric Calabi-Yau n-fold is represented as a T n fibration over an 
n-dimensional polyhedral cone C. The boundary of C consists of (n — l)-fans. On 
each (n — l)-fan a cycle v in the toric fiber shrinks. In other words, the fan is 
the fixed submanifold of the U(l) isometry generated by the vector v. For each 
(n — l)-fan, there is a vector v representing the shrinking cycle, and the toric 
data is given as a set of such vectors. 

In order to extract the toric data of the moduli space from the information 
of a gauge theory, it is convenient to translate the system into a gauged linear 
sigma model (GLSM). This is achieved by solving the F-term conditions with the 
help of perfect matchings. 

A perfect matching is a number assignment /// to links in a tiling which 
satisfies the following conditions. 

• fij = or 1 for any link I. 

• Among links ending on a vertex a, only one has non- vanishing /j. 
The following equation follows from these two conditions. 

(7»,A0=E/r = 1 Va ' ( 21 ) 



7 



fi* is the tangential flow associated with the number assignment fij, and the 
product (*, *) is the intersection of a normal flow and a tangential flow, which is 
defined by 

(f,9*)=Y,fi9i- (22) 
Figure H] shows the four perfect matchings of the AB JM tiling. 




Mi /x 2 /x 3 



Figure 4: The four perfect matchings of the ABJM tiling. 



The F-term conditions require all the terms C 7a in the superpotential are the 
same. We can solve this condition by [55] 

*z = II P» ( 23 ) 

where is a GLSM field defined for each perfect matching /i, and the summation 
is taken over all the perfect matchings with \ii = 1. Substituting this into the 
definition of la and using (!2il . we can show that C-y a is the product of all the 
GLSM fields regardless of the index a, and the F-term conditions are therefore 
satisfied by (1231) . Because this expression is redundant, we need to extend the 
gauge symmetry G acting on by adding £7(1) rotations of GLSM fields which 
keep invariant. Let G' be this extended gauge symmetry. If the number of 
the perfect matchings is n pm and the space spanned by the GLSM fields is C npm , 
the moduli space of the gauge theory is the coset C npm /G'. 

In order to obtain the toric data, we need to find U(l) symmetries which 
have non-trivial fixed submanifolds. It is easy to show that in the moduli space 
defined as the coset C Upin /G', such a submanifold is given as the image of the 
fixed plane p M = of 11(1)^ symmetry by the homomorphism C rapm — > C Upm /G', 
where 11(1)^ is the symmetry which rotate only one GLSM field p M with charge 
1. Thus, the components of the killing vector are given as the 11(1)^ charges 
of the four toric coordinates. 

As we mentioned above, we can use the four gauge invariant operators in fl20l) 
as coordinates in the Calabi-Yau 4-fold, and then the four components of are 
C/(1) M charges of these operators. By substituting fl23|) into (@J we rewrite an 
operator O q in terms of GLSM fields as 

o,=nn#=iip? ,M * > > ( 24 ) 

I H3l M 



s 



and thus, the U(l)^ charge [O q ]^ of the operator O q is given by 



(25) 







f («,**•) \ 






(Am*) 












V i / 



Applying this formula to the four gauge invariant operators, we obtain the fol- 
lowing components of the killing vectors 



(26) 



The last components of these vectors are always 1, and this guarantees that the 
toric manifold is Calabi-Yau. With this formula, we can easily obtain toric data 
from a given IIA tiling. 

As a simple example, let us consider the ABJM tiling in Figure HJ If we use 
the four perfect matchings in Figure IU and flows s, ct, (3, and j a in Figured] and 
[31 we obtain the following four killing vectors. 



(27) 



By neglecting the fourth components of these vectors and plotting corresponding 
points in the three-dimensional lattice, we obtain the toric diagram of the moduli 
space C 4 /Z fc . (Figure E]) 
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Figure 5: The toric diagram of the orbifold C 4 /Zk. 



6 Relation to crystals 

A brane tiling describing a four-dimensional gauge theory can be regarded as a 
brane systems consisting of D5-branes and NS5-branes, which is T-dual to D3- 
branes probing a toric Calabi-Yau 3-fold, and the rules of reading off the gauge 
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Figure 6: 2-cycles used to extract toric data from a crystal are shown. These are 
expected to represent M2-branes corresponding to gauge invariant operators! 



theory from the tiling have natural interpretation in terms of this brane system. 
For example, faces in a brane tiling represent D5-branes in the brane system, and 
the U (N) factors in the gauge group are identified with the gauge groups realized 
on the D5-branes. 

Brane crystals [29| 1301 |3T] are analogues of brane tilings for M2-branes prob- 
ing four- dimensional toric CY cones. By T-duality transformation in M-theory, a 
system of M2-branes probing a four-dimensional toric Calabi-Yau cone is trans- 
formed into a brane system consisting of M5-branes. Brane crystals are bipartite 
graphs in T 3 representing the structure of the M5-brane systems [29 J. 

Contrary to the case of brane tilings for four- dimensional gauge theories, we 
can obtain much less information from this brane system. This is because we 
have only little knowledge about theories realized on M5-brane systems. The 
purpose of this section is to obtain some information about the relation between 
Chern-Simons theories and brane crystals by using the results obtained in the 
previous sections. 

The method to obtain toric data from crystals has been already known[30j. 
Actually, we may define brane crystals as bipartite graphs in T 3 which give toric 
data of Calabi-Yau 4-folds in a similar way as brane tilings. The method to 
obtain the toric diagram from a brane crystal is as follows: First, instead of the 
a and (3 cycles in brane tilings, we define three closed 2-cycles A (1-3 plane), B 
(2-3 plane), and S (1-2 plane). (See Figure [6j) We assume these do not include 
vertices of the crystal on them. We also define a cycle G a for each vertex a, 
which is a sphere enclosing the vertex a. (See Figure [6j) Because the crystal is 
bipartite as well as tilings, we can define perfect matchings on it. If we denote 
the intersection number of a 2-cycle C and a perfect matching fi by [C, /i] , the 
vectors forming the toric diagram are given by 



( [A/i] \ 
[S, fi] 

V [G a ,i4 J 



(2* 



By definition, the last component [G a ,fj] is always 1, and the Calabi-Yau condi- 
tion is satisfied. 
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Figure 7: A crystal for ABJM model with k = 2 is shown. This figure includes 
two fundamental regions. 

When a IIA tiling is given, it is easy to construct a brane crystal which gives 
the same toric data by the formula ( 1281) as the data obtained by ( 1261) from the 
IIA tiling. Let ( 

XaiDa) be the coordinates of the vertex a in the IIA tiling. We 
put the corresponding vertex in the crystal at the point (x a ,y a ,0) in the three- 
dimensional torus. We make the three-dimensional graph by connecting these 
vertices in the same way as the tiling. Namely, if vertices a and b in the tiling are 
connected by a link, we connect the corresponding points in T 3 by a link, too. 
There are infinitely many ways of connecting two vertices in T 3 with a link with 
different winding numbers. We fix this ambiguity by requiring the following two 
conditions. 

• The crystal reduces to the original tiling by the projection along the vertical 
axis. 

• The vertical winding number of link / is s/. 

In other words, we interpret the integers si assigned to links as the gradient of 
links in the three-dimensional space. See Figure [7J for an example of the ABJM 
model with k = 2. As the result, we obtain a three-dimensional bipartite graph 
with the same number of vertices and links as the original brane tiling (Figure [7]). 

Let us confirm that the crystal constructed in this way correctly reproduces 
the toric data ( |26i) obtained in the previous section. First of all, the two- 
dimensional and three-dimensional graphs are differ only by their embeddings 
to the tori. The former is embedded in T 2 and the latter is in T 3 . Therefore, the 
three-dimensional graph has the same perfect matchings as the two-dimensional 
one. 

Let us first consider the first two components (a, ^i*) and (f3, fj,*) in (1261) . We 
define the three-dimensional lift of a and (3 cycles by 

A = a®S\ B = [3®S\ (29) 
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where S 1 here is the cycle along the vertical direction. We use these 2-cycles as 
A and B in the formula (1281) . Then it is obvious that the first two components 
of ( l26l) and those in (|28l) are the same. 

(a,n*) = [A,fi], (J3,n*) = [B,n]. (30) 

For the last components, we define 2-cycle 7 a ® S 1 for each vertex a. These 
are homologous to G a defined above, and the relation 

(7«,M*> = [G a ,d (31) 

holds. (Actually, this is by definition always 1.) 

Finally, let us consider the third component in (1261) . which is given as the 
U{l)n charge of the baryonic operator B. In the brane tiling, the baryonic op- 
erator is expressed in different way from the other mesonic operators. Mesonic 
operators are associated with conserved flows on the tiling, while the flow s cor- 
responding to the baryonic operator B is not conserved. In the crystal, however, 
the third component is also given as the intersection of closed 2-cycle and perfect 
matchings. As we mentioned above, the link / in the crystal has the vertical 
winding sj, and it intersects sj times with the 2-cycle S. Therefore, the intersec- 
tion (s, /j,*) can be rewritten as the intersection of the closed 2-cycle S and the 
perfect matching /i. 

(«,***> = tM- (32) 

Thus, the third components of (1261) and ([28]) are the same. 

Now we have confirmed that the crystal constructed above correctly repro- 
duces the toric data of the moduli space of the Chern-Simons theory described 
by the tiling. At the same time, we have established the correspondence between 
the gauge invariant operators in (I2UI) and closed 2-cycles in the crystal. Because 
the set of operators in (lT3"j) and (IT31) generates arbitrary gauge invariant opera- 
tors, we have established the complete map between gauge invariant operators 
including both mesonic and baryonic ones and closed 2-cycle in the crystal, which 
are interpreted as closed M2-branes[29[ 130] . An interesting feature of this cor- 
respondence is that even though in the Chern-Simons theory baryonic operators 
and mesonic operators have different structure, the brane crystal describes these 
in the parallel way. 

As another support to our prescription, we can show that the level kj are 
naturally obtained from the brane system described by the crystal. In order to 
read off the Chern-Simons theory from a crystal, we need to project the crystal 
along the vertical direction, and go back to the tiling. From the viewpoint of brane 
system, we can interpret this projection as the compactification of M-theory to 
type IIA string theory. Then, links and faces in the tiling are interpreted as 
a network of NS5-brane and D4-branes ending on the NS5-brane, respectively. 
Gauge groups are realized on the D4-branes, and the Chern-Simons terms are 
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induced from the following boundary term in the D4-brane action. 

S =— I A A dA A d(f), (33) 

2lX JdT>4 

where A is the gauge field on the D4-brane and 4> is the compact scalar field on 
the NS5-brane corresponding to the X 11 coordinate in the M-theory picture. If 
the scalar field <fi has non-trivial profile along the boundary of the D4-brane, this 
induces the Chern-Simons coupling in the three-dimensional gauge theory, and 
the level is given by 

ki=fd<f>, (34) 



where the integration is taken over the boundary of the face i. If we identify d<p as 
the gradient sj along links, fl34"|) is nothing but the relation ([2D, or, equivalently, 
(El). 



7 Conclusions 

In this paper we investigated the relation between brane tilings describing M = 2 
Chern-Simons theories and the toric data of their moduli spaces. We gave a simple 
procedure to read off the toric data of the moduli space from the brane tiling. 
In order to obtain the toric data, we should first represent the Chern-Simons 
levels as a flow s on the tiling, and the vectors forming the toric diagram are 
obtained as the intersection of the perfect matchings /x* and the flows (a, (3, s). 

IIA brane tilings, which are regarded as brane systems consisting of D4-branes 
and NS5-branes, can be regarded as the projection of the crystals, which describe 
M5-brane systems. We showed that we can lift a IIA tiling to the corresponding 
crystal by using the flow s* as the gradient of links. We found that gauge invariant 
operators, which include both mesonic and baryonic ones, are represented in the 
crystal as closed 2-cycles. 

We emphasize that although our prescription always gives a crystal for a 
given IIA tiling, it is not always possible to give a tiling which reproduce a 
given crystal. Our prescription does not guarantee the existence of a Chern- 
Simons theory which reproduces a given toric Calabi-Yau 4-fold as its moduli 
space. There may not exist corresponding Chern-Simons theories for a class of 
manifolds. Contrary, there are crystals which gives more than two tilings by 
the projection along different directions. This may suggest the duality among 
Chern-Simons theories. We wish to come back to these issues in near future. 
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